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Abstract 

We review developments made since 1959 in the search for a closed 
form for the susceptibility of the Ising model. The expressions for the 
form factors in terms of the nome q and the modulus k are compared and 
contrasted. The A generalized correlations C[M,N;\) are defined and 
explicitly computed in terms of theta functions for M = A'^ = 0, 1. 

1 Introduction 

There are three important thermodynamic properties of any magnetic system 
in zero magnetic field: the partition function from which free energy and the 
specific heat are obtained; the magnetization; and the magnetic susceptibility. 
For the two dimensional Ising model in zero field defined by 

£a = —'^^{E^aj^kcrj+i,k+E^(yj,k(yj,k+i} (1) 

with aj^k = ±1 the free energy was first computed by Onsager [T] in 1944 and the 
spontaneous magnetization was announced by Onsager in 1948 [5] and proven 
by Yang [3] in 1952. To this day a closed form for the magnetic susceptibility has 
never been found. Wc will here trace the saga of the quest for this susceptibility. 

If we could solve the Ising model in the presence of a magnetic field H which 
interacts with the total spin of the system as 



£ — £q — Hy^aj^k 



(2) 



then the magnetic susceptibihty would be computed as 



^ 'Jip (3) 



where the magnetization is 



M{H) - Jtt, i: ^o^e-'/'^^ (4) 

with the partition function defined by 

Z{H) = e-^/''^^ (5) 

However, because the Ising model has only been solved for = we are forced 
to restrict our attention to x(0) which from ([2])- ([51) is given in terms of the two 
point correlation functions as 

keT ■xiO)= Y.{{ao,o<JM,N)-MiOf} (6) 

M.N 

where M(0) is the spontaneous magnetization of the system which is zero for 
T>Tc and for T < 

M(0) = (1 - fc2)i/8 (7) 

where 

k ^ (sinh2X^'sinh2if'*)^i (8) 
with if"''' = E^'^^/kBT and is defined by 

fc = 1. (9) 

The first exact result for the susceptibility was given in 1959 by Fisher [3] 
who used results of Kaufmann and Onsager [5] to argue that as T — > Tc the 
susceptibility diverges as \T — Td"^^'*- The saga may be said to begin with the 
concluding remark of this paper: 

In conclusion we note that the relatively simple results (1) and (3) suggest 
strongly that there is a closed expression for the susceptibility in terms of elliptic 
integrals. It is to be hoped that such a formula will be discovered,- ■ ■. 

One year later Syozi and Naya[6], on the basis of short series expansions, 
proposed such a formula for T > which does not involve elliptic integrals 

(1 - sinh^ 2/s:^ sinh^ 2/s:'*)i/4 
^ '^^>~ cosh cosh 2if - sinh 2K^ - sinh 2K^' ' 

However, when higher order terms were computed this conjecture was shown 
not to be exact. 

To this day the "closed expression" for the susceptibility hoped for in [4^ has 
not been found. 



2 



2 Form factor expansion and the A extension 



To proceed further a systematic understanding of the two point correlation 
function is required. For short distances the correlations are well represented by 
determinants OE] whose size grows with the separation of the spins. However, 
in order to execute the sum over all separations required by © an alternative 
form of the correlations which is efficient for large distances is needed. The 
study of this alternative form was initiated in 1966 by Wu [5] who discovered 
that for the row correlation (cto.oCo.w) that when 
TV - IT-Tcl > 1 for T <Ta 

{ao,o<Jo.N) = (1 - ty^' ■ {1 + /S + • • •} (11) 

with t ^ P with k given by dH). For T > we define k = sinh2i^" sinh2is:'' 
and find 

(ao,o^o,iv> = (1 - t)'^' ■ {f!,% + ■■■} (12) 

with t = k'^. In (|lll) and ([T^ /^'^ is an n fold integral which exponentially 
decays for large N . The results ()TT|) and (IT^ are the leading terms in what 
has become known as the form factor representation of the correlations which 
in general for T < Tc is 

oo 

(ao,o^M,^) - (1 - t)^" ■ {1 + E /Si (13) 



and for T > 



OO 



KoCXM,^) = {l-tf"-Y.fMT^ (14) 

n=0 

where /|"''^ is an n dimensional integral. For general M, N these were 
computed in 1976 by Wu, McCoy, Tracy and Barouch [9l and related forms 
are given in |101 [TT| . However, for the diagonal correlations an alternative and 
simpler form is available which was announced in |12| and proven in [13,. For 
the diagonal form factor for T < Tc 



^niN+n) ^ ^ ^ ( (l-fa2,-)(x-/-l) \ 

fcl 7=1 V (1 - *^2j-l ) (2:27-1 - 1) / 

2 

2 



(15) 



l<j<ri l<fe<ri ^ ^/t 1 / l<j<A;<?i 



for T > 



f(2n+l)/,N 
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^{n+l/2)N+n{n+l) rl 2n+l n+1 
■ V / ■ ^ f — 1 ■ _ 1 



k=l J=l 

1 



X2k 



n-2,[(i-te,,)(.,7-i)]^/^ n n i^^^ 

j = \ l<j<n+l \<k<n ^ ^ 

[xij-i - X2k-\f (a;2j-a;2fe)^ (16) 

\<i<k<n+\ l<j<k<n 

In particular 

where -F(a, 6; c; is the hypergeometric function. 

It is often useful and instructive to extend the form factor expansions ([1 



and by weighting /|"''^ by A" and thus we define "A generalized correla- 
tions" 

oo 

C_(M, N; A) - (1 - . {1 + ^ A2"/(^n) } (18) 

n=l 

and for T > 

OO 

C+(Af, AT; A) = (1 - t) . ^ a2"+V(?"+i) (19) 

n=0 

This A extension was first introduced in 1977 by McCoy, Tracy and Wu [14] in 
the context of the scaling limit. 

3 Leading divergence as T — > Tc 

These form factor expansions may now be used in ([6]) where the sums over 
M, N are easily executed under the integral signs to produce a corresponding 
expansion of the susceptibility fS] which we write for T < as 

oo 

ksT-xiO) = (l-O'/'-Ex^'"^ (20) 

Tl=l 

and for T > Tc 

OO 

ksT -xiO) = (1 - . ^ (21) 

where 

oo oo 

X'"'- E E (22) 

_A/=-oo N=-oo 
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For n = 1,2 the x^"^ are explicitly evaluated [5]. For the isotropic lattice we 
have 

^ 3^(l-fc)(l-A:2) 

where K and i? are the complete elliptic integrals of the first and second kind. 
The form (jTU]) of Syozi and Naya 6 with = is seen to be the first term 
in (mi). It is quite clear that unlike the 1959 argument of 5] the behavior of the 
susceptibility as T —> Tc will be different depending on whether T approaches 
Tc from above or below. This dramatic difference was first seen in 1973 in [T5] 
where it is shown that for T — s> Tc± 

kB-x{0) - Co±-\T-T,\-'^^ (25) 

where 

Co± = 2-1/2. coth2i^^ coth2ii:^-[i^^' coth2ii'^" +ii:^^coth2i4:^]-^/'^-/± (26) 

and I± have been numerically evaluated to 52 digits in |16) : 

/+ = 1.000815260440212647119476363047210236937534925597789 • ■ • 

J_ = • 1.000960328725262189480934955172097320572505951770117 • • • 
127r 



4 The singularities of Nickel 

The next advance in the understanding of the analytic structure of the suscep- 
tibility came in 1996 when Guttmann and Enting |17) . using resummed high 
temperature series expansions of the anisotropic Ising model, argued that the 
susceptibility cannot satisfy a finite order differential equation and raised the 
question of the occurrence of a natural boundary. This natural boundary ar- 
gument was made very concrete for the isotropic case in 1999 and 2000 by 
Nickel [18] and [19] who analyzed the singularities of the n fold integrals x''"''- 
These integrals, of course, have singularities at T = Tc where the individual 
correlation functions {(Jofi<^M,N) have singularities. However, Nickel made the 
remarkable discovery that the integrals, for x^"\ contain many more singulari- 
ties. In particular he found that, for the isotropic lattice, x'""* has singularities 
in the complex temperature variable s = smh2E/kBT at 

s = Sj,k = e'''-' (27) 

where 

2cos(6'j^fc) = cos(27rfc/n) -I- cos(27rj7n) (28) 
For n odd (T > T^) the behavior of x*-"-* near the singularity is 

^(2„+i) _ ,2«("+i)-i.ine (29) 
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with 

e = 1 - s/sj^k (30) 

and for even n (T < Tc) 

The discovery of these singularities demonstrates that the magnetic sus- 
ceptibihty is a far more comphcated object than either the free energy or the 
spontaneous magnetization and that the hope expressed in [4] of a closed form 
in terms of a few elliptic integrals is far too simple. 



5 The theta function expressions of Orrick, Nickel, 
Guttmann and Perk 

In the following year a major advance was made by Orrick, Nickel, Guttmann 
and Perk [1^1 who studied both the form factors and the susceptibility by means 
of generating on the computer series of over 300 terms. From these series they 
then made several remarkable conjectures for the form factors. 
To present these conjectures we define theta functions as 

CXj 

ei{u,q) = 2 ^ (_i)"g(»+i/2)^ . sin[(2n+ (32) 

n=0 

oo 

02{u,q) = 2 J2 9^"+'/'^' • cos[(2n + l)u] (33) 

n=0 

oo 

03{u,q) = 1+2^ g"' •COS271U (34) 

n=l 
oo 

6l4(u,q) = 1 +2^(-l)"g"' •cos2nM (35) 

n=l 

and for u = we use the short hand 

02^ 0^(0, q), 03(0, q). 0i^0i{{),q) (36) 

The quantity q is the nome of the elliptic functions and is related to the 
modulus k by the relation 



l + q' 



(37) 



In terms of these theta functions, conjectures for form factors are given in 
sec. 5.2 of [1^ by defining an operator $o which converts a power series in z to 
a power series in q as 

oo oo 
n=0 n=0 
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Conjectures are then given for /o'o\ fi"i, f20 and /j"-*. In particular 
we note 

2- • (1 - k^y^^ • = . ^^YT^) ^''^ 

and 

where 

(l,A;-^/2)= 1 for r<rc (n even) 

foj. T>rc (nodd). (41) 



6 Linear differential equations 

A second approach to the form factors and susceptibihty was initiated in 2004 in 
[20] and subsequently greatly developed in [HI [11 [23l [Ml [25l [26] . These studies 
are similar to '16' in that they expand the form factors and susceptibility in long 
series. However, instead of the nome q the expansion is in the (modular) variable 
t. The goal of these studies is to characterize the n particle contributions x'"^ (t) 
to the susceptibility in terms of finding a Fuchsian linear ordinary differential 
equation satisfied by x^^Ht)- Such a linear differential equation always exists 
for an n-fold integral with an algebraic integrand in some well-suited choice of 
integration variables and in the parameter t. However the order and the degree 
of the equation rapidly become large for increasing n and it may take series of 
many thousands of terms to find the differential equation. Such a study can 
only be done by computer. 

There are several features of these differential equations to be noted. In 
particular the operator which annihilates x^""* factorizes and furthermore the 
operator has a direct sum decomposition such that x^""^-'-' for j = 1, • ■ • [n/2] 
are "contained" in x'""*- 



7 Diagonal form factors 

With the observation of factorization, direct sum decomposition and Nickel 
singularities of the n particle contributions to the bulk susceptibility x^"''j it- 
has become clear that the susceptibility is far more complicated than what was 
envisaged by Fisher ^ in 1959. Because of this complexity the question was 
asked if there could be a simpler object to study which would yet be able to give 
insight into the structures which had been observed. Several such "simplified" 
objects have been studied [27] which consist of more or less forcibly modifying 
parts of the integrals for the x*'"'- However, there is one "simplified" model 
which commands interest in its own right. This is the "diagonal susceptibility" 
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which is defined [55] by restricting the sum in © to the correlation of spins on 
the diagonal 

oo 

kBT-Xd= ^ {Koajv,Jv)-Af2(0)}. (42) 

JV=-oo 

In statistical language this diagonal susceptibility is the susceptibility for 
a magnetic field interacting only with the spins on one diagonal. In magnetic 
language this is the p — value of the groundstate structure function 

oo 

S%p)= ^ e'^'HK'^J)-^!} (43) 

j = -oo 

of the transverse Ising model 

oo 

Hti = E {^."^J+i+^'^D (44) 

j = -oo 

These interpretations give the diagonal susceptibility a physical interpreta- 
tion which the other "simplified" models do not have [27]. Furthermore much 
more analytic information is available for the diagonal Ising correlations than 
for correlations off the diagonal. Firstly it is known from the work of Jimbo 
and Miwa [29^ that the diagonal correlations are characterized by the solutions 
of a particular sigma form of Painleve VI equation and secondly the integral 
representation of the diagonal form factors (IT5|) and (|16l) is more tractable than 
the representation of the general off diagonal correlations. 

The diagonal form factors have been extensively studied in [T^ by means 
of processing the differential equations obtained from long series expansions by 
use of Maple. Diagonal form factors /^]v f*^'" " ^'^ large as 9 and N as large 
as 4 have been studied and many examples are given in |12| where they have 
all been reduced to expressions in the elliptic integrals E and K. A few such 
examples are as follows: 

For n = 1 (when the hypergeometric function of (fTTl) is reduced to the basis 
of E and K by use of the contiguous relations) 

fl,]^ = {2M ■ K (45) 

^1/2 /(I) ^ (2/7r) -{K^E} (46) 

3tf^]^ = {2/n)-{{t + 2)K -2{t + l)E} (47) 

15 ^3/2 f^^^ = (2/7r) • {{4t^ + M + 8)K - {8t^ + 7t + 8) E} (48) 
105^2/1,4 = {2/Tr) ■ {{241^ + + mt + 48)K 

- {m^ + + m + 48) i;} ; (49) 

2/(^J = {2/nr -KiK-E) (50) 



for n 
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2 f['^l = 1 - (2/7r)^ ■K-{{t-2)K + 3E} 
(2/7r)2 . {6^2 _ lit + 2)K'^ + (15t -~4)KE + 2{t + 1) E^] 



3150 tVi'i 



135t2 - (2/7r)2 . {(137i^ - 242^^ 



52t- 



319r + 122t + 16) KE + 4(t + l)(2t^ + 13t + 2)E^) 



(51) 

(52) 
(53) 
(54) 



6300 

~{2/iif ■ {{m^ + 6440t'' - 1119t^ + 2552t2 + 464t + 128) K"^ 
-{12U^ + 576*"^ - 14519t^ + bAU"^ + lG56t + 256) KE 
+ (1 + i)(16i'' + 58i^ + 333i^ + 58t + 16) E"^}; 

for n = 3 

6 /^'o^ = (2/7r) • K - (2/Tif ■ K^-{{t-2)K + iE} 



(55) 



(56) 

6^1/ Via = A{2/ti) ■ {K - E) - {2/iTf ■ K ■ {{2t-i)K'^ +QKE -iE"^} 

(57) 
(58) 



18i/^^2^ ^ 7 (2/7r) ■ {(t + 2)i^ - 2(t + 1)£;} 

iii + 2)ii:2i; 



(59) 



- (8i2 + 7t + 8)i£;} 

if 3 



-(2/7r)3 • {3(t2 - 2) if^ - 3 {2f 
-36 (t- 1) A'£;2 - 24£;3} 
270i^/2j^3) ^ 30 (2/7r){(4i2 + 3i + 8)if 
-(2/7r)3 • {(72t^ - 158<^ + 189i^ - 156t + i , 
-6(24t^ - 108^"^ + 2Qt^ - 6t + 4) K'^E 
-3(232*3 -111*2-180* -8)KE^ 
-4 (t + 1)(2<2 + I03i + 2) t E^}; 

for n = 4 

24 /(^J = 4(2/7:)' -KiK-E) 
-i2/nf ■ /s:2{(2< - 3) + 6KE - 3E'} 
24f['^l = 9 - (2/7r)2 . 10/s:{(t-2)/s: +3£:} 
+ (2/7r)* ■ /s:2{(t2 -Qt + 6)K' + 10 (t - 2) ifi; 
72t/i^2 = 72t 

-(2/7r)2 . 16 • {(6*2 - lit + 2)K' + (15< -4)KE + 2{t + 1) £;2} 
+ (2/7r)4 • {24t3 - 98f2 + im - 36) X'' + 2 [74^ - 157t + 66) K^E 
+3 (71i - 60) K'^E' + 12 (i + 9) KE^ - 24E^}. (63) 

These examples are sufficient to illustrate the following phenomena which 
hold for all examples considered in [12] and which are certainly true in general: 



IbE'} 



(60) 

(61) 
(62) 
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/i.'"^ = ll^i,n-9%%{t) (64) 

n 

= T.4n-9'^T\t) (65) 



J=0 

T 

f(2n+l) - ^ -+ - + 
3=0 

where c^„ are constants independent of t and for even j are of the form 



9So\t) = E^^;oW-^'""''^'' (66) 

j=o 

n 

9?,i\t) = E^^a(^)-^'""''^' (67) 

2n 

9N$it) = t-""^' E Pj'nM*) ■ ^ ^ 2 (68) 

and for odd j 

n 

9^T'\i) - E ^^^.o(^) • (69) 

71+1 



(2n+l) 



(t) = t-'/' ■ E P^n;lit) ■ (70) 



(2n+l), 



2ri+l 



ffiv.^ '(i) = t-''/' ■ E ^^.^(0 • K'''^'-' E' for ^ > 2 (71) 

where Ppn-mi^) ^'^'^ polynomials. 

The decompositions (|64)) and (|65|) represent a direct sum decomposition of 
the form factors [28j . The fmictions 5]v V individually are annihilated by Fuch- 
sian operators which are equivalent to the j + 1 symmetric power of the second 
order operator assosciated with the complete elliptic integral E (or equivalently 
K). 

We also observe the relation between f^i{t) and fj^"f^'^\t) 

(2/^) . K ■ f[f\t) = (2n + 1) . /^'o"+''(t) (72) 
{2/n) . . K ■ f[',^+'\t) = 2 (n + 1) . f!,f+'\t) (73) 
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8 Nome ^-representation versus modulus /^-representation 

We will need the following identities which relate functions of the nome q = e""^ 
where r = iK{k') / K(k) with functions of the modulus k 

k= |, k'^{i-eY'^^ |, (74) 



n 

which we will use as 



-K^e,. and -^-^^ (75) 



We will also use 

dK _ E-k''^K dE _ E-K 

Ik " AJP ' Ik ~ k ^ 

8.1 fif 

We first write for j = 2n using (TMl) as 

(2„) _ 1 / 2^"z^".(l-z^) N 
/o,o - ^-^o^ (i + ^2)2„+i J (78) 

Thus, by use of the elementary expansion 



(1 + Z2)2n+1 (2„), 



and the definition (|38p of the operator we find that in terms of the nome q 



J=0 



n-1 



(-1)"4" A, TT r ^ 



6'4(2n)! ' J--'-''' do 

(-1)"4" d 21 . 

To convert this to an expression in terms of the modulus k we first use (|74p 
to write 

■k' -K (81) 



9 2 

^4 = 

TT 



and thus using ([75]) 

' dq ^ TT^ dk \t: 



q^0l ^ l^.kk'^-K^-4r(-k'K] (82) 
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which using ([77| reduces to 

dq TT'^ TT 

Using ((8T|) on the right hand side we find 

d „ 2 



9—^4 = — ■ k' ■ ■ -{E - K}. (83) 



204 • <Z:f ^4 = ■^■0l-K-{E~K} (84) 
and thus 

^-K-iE-K}. (85) 

To evaluate /q^q we use (|85|) in ((80|) with n 1 to obtain 

/o'o = 4 • if • {i^ - i?} (86) 

TT 

which is in agreement with (j50p . For arbitrary n the form factor /g^J*^ is obtained 
from dHO]) by repeated use of dH]) and ([771). 

e D /.(2n+l) 

o-^ Jo,o 

To study /^Jo"^^^ we first use ^ to write dMl) as 

- (l + z2)2n+2 j (87) 



and then, using the elementary expansion 

22»+l . ^2n+l . , ^2) 
(1 + z2)2"+2 

= 2 i— • J2 (2J + 1) • (-IF • • n [(2.? + 1)' - (2"^ + 1)' 

and the definition (1551) of the operator $o we find 



(88) 



j^(2n+l) 

/o,o 



73 



n-l 



'3 2( 1)" J|[4,jL_(2„^+l)2]^(2^- + l)(-l).g(2.-+i)V|89) 



Thus, if we write 



2^(2j + l).(-l)^.g(2, + l)V4 ^ ^ g^g^g^ (gg) 



9w 
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where in the last hne we have used a well known identity, we find the result 

(-1)" 'fr^ ^ d 



We may now use (rMl) - ([77|) to reduce (l^lT) from a function of q to a function 
of yfc. 

For n = we use (|75p to find 

/^,o^ ^ el= ^-K (92) 

TT 

which agrees with (j45|) . 

For n > 3 we need an expression analogous to ((55)) for the product 62O394. 
From (Ull), (123 



6*20202 = ^fc' . (2X/7r)3 (93) 

and thus 

5^0203^1 = - A.fcfc'2i^2.i.{fcA:'(2X/7r)3} 

= A;-kk' ■ (2K/tt)^ ■ K ■ {(fc2 ~2)K + iE} 

= ^■^2^3^4-^-{(^'~2)X +3S}. (94) 
Therefore we obtain 

q^e^e^e^ = ^02^3^4 • k ■ {{e -2)k + 3^} (95) 

which when used in (j9ip with n — 1 gives 

- ^ • {(2A) • ^ - (2/7r)3 if • [(fc2 ~2)K+ iE]} (96) 
which is in agreement with (j56p . 

8.3 

The equalities ([72]) and (|73l) which express in terms of /q'^'"' follow imme- 
diately from ([39]) and (00]) by use of ^ and ([25]). 



9 The A generalized correlations 

For iV = 0, 1 the diagonal A generalized correlations defined by (|18l) and (|19p 
may be obtained by using the expressions for (1 — t^^'^f^j^^ in terms of the 
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operator $o as given by ([5^ and (HUl) . In this form the sums over n are easily 
done as geometric series and the linear operator is then used to convert the 
series in z to series in the nome q which can then be expressed in terms of 9 
functions as was done in the previous section. Then, setting 

A = cosM (97) 

we obtain the foUowing resuUs 

C.(0.0;A, . „ 



= (101) 



siTi(>i)fli(0;9)93(0;?) 

where prime indicates the derivative with respect to u. The result (|98l) was first 
reported in jl2]. The resuhs (|99p - (|10ip have recently been given in ^30^. For u a 
rational multiple of tt these diagonal generalized correlations reduce to algebraic 
functions of the modulus k. Several examples for C_(0,0;A), C_(1,1;A) and 
C_(2, 2; A) are given in [12]. 



10 Diagonal susceptibility 

We may now explicitly obtain the diagonal susceptibility by using the form 
factor expansion p^ - (IT6|) in the definition p2)) and evaluate the sum on N as 
a geometric series. We obtain for T < Tc 

oo 

kT-x,-{t) = {l-ty/*-Y.xT\t) (102) 



with 



^'^^ w^-X ■■■Jo n/^-i-t".....x.„ 

^TT f X23-l{l-X2j)il~tX2j) \ ^'"^ 

7=1 V^2j(l - a;2i-i)(l -1x2^-1)) 

X Jl n (1 - ^ 2^2^-1 2^2*;)^^ 
\<i<n l<k<n 

X n {X2j-1 - X2k-l)'^ {X2j - X2k)'^ (103) 
l<j<fe<n 
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and for T > 

OC 

kT ■ xa+it) = (1 - tf'^ ■ J2 xT^'\t) (104) 



n=0 



with 



1 +in+l/2^^...^2^^^ " , .1/2 



n+l 

n -2;2i-i)(l -<X2j-i) •a;2j-i 



X 

^ n n (1 -^a;2j-ia;2fe)-' (105) 

l<j"<n+l l<k<n 

l<j</£<n+l l<j<k<n 

This diagonal susceptibihty has been extensively studied in [55] • 
The integrals for xll\t) and x|i'^''(t) are explicitly evaluated as 



~(i) 

and 



tiit) = Y~W ^^^^^ 



Fuchsian equations have been obtained for xf\t), x'd\t) and x''d\'^)- 

(3) 

From these equations we find that jQ (t) has a direct sum decomposition 
into the sum of three terms. One term is just Xd\t) given by (|106|) : the 
second is 

and the three solutions to the differential equation for the third term are two 
Meijer G functions and 

(l + 2fc)(fc + 2) 



(l-fc)(l + fc + /c2) 



{^^(1/6,1/3;1;Q)^ 



+ ^ • F(l/6, 1/3; 1; Q)F{7/6, 4/3; 2; Q)} (109) 

where 

^ 4(fc2 + fc+l)2 ^^^^> 
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Furthermore the Xd (t) have singularities on \t\ = 1 which are the analog of 
the Nickel singularities for the bulk susceptibility. For T < Tc the singularities 
in x^^^HO at = 1 and are of the form e^" ~^lne and for T > Tc the 
singularities in ;^(2ri+i) ^^j-g ^n+1/2 _ ^ and are of the form e("+i) 

11 Natural boundary 

The most intriguing feature in both the bulk and the diagonal susceptibility are 
the singularities on the unit circle of the modular variable k = 1. As n increases 
the number of these singularities increases and becomes dense as n — 00. 
Therefore, unless a massive cancellation occurs the susceptibility will have a 
natural boundary on the circle |fc| = 1. Recently further arguments in favor of 
such a natural boundary were given in ^23j. In terms of the nome q the circle 
|/c| = 1 corresponds to the curve in Fig. 1. 




Figure 1 : The curve in the plane of the nome q of the unit circle | /c | = 1 



12 Conclusion 

Wc have seen that since 1959 a great deal of progress has been made in under- 
standing the susceptibility of the Ising model and that the analytic structure is 
vastly more complicated than was envisaged 50 years ago in [4]. In particular 
the existence of a natural boundary is a completely new phenomenon which has 
never before appeared in the study of critical behavior. The connections with 
elliptic modular functions are profound and extensive and much of the struc- 
ture still remains to be discovered. It is quite remarkable that in 50 years the 
problem has not been solved. 
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